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Abstract. We combine ab initio density functional theory (DFT) structural studies 
with DFT-based nonequilibrium Green function calculations to investigate how the 
presence of non-hexagonal rings affects electronic transport in graphitic structures. 
We find that infinite monolayers, finite-width nanoribbons and nanotubes formed of 
5-8 haeckelite with only 5- and 8-membered rings are generally more conductive than 
their graphene-based counterparts. Presence of haeckelite defect lines in the perfect 
graphitic structure, a model of grain boundaries in CVD-grown graphene, increases 
the electronic conductivity and renders it highly anisotropic. 
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1. Introduction 

Graphene is a unique 2D material that combines extraordinarily high electrical and 
thermal conductivity [1, 2] with mechanical strength, flexibility, thermal and chemical 
stability. Interest in this system increased significantly after a successful mechanical 
exfoliation using a Scotch tape has been reported [3] that yielded large, defect-free 
samples. As a scalable alternative to the “Scotch tape” exfoliation technique, chemical 
vapor deposition (CVD) is commonly being used now to form graphene monolayers on 
metal substrates including Cu [4, 5, 6, 7]. The quality of CVD-grown films is inferior to 
those produced by exfoliation, since gas phase deposition leads to simultaneous growth of 
graphene flakes that eventually interconnect, forming grain boundaries with a defective, 
haeckelite-like structure [8, 9, 10, 11, 12, 13] consisting of non-hexagonal carbon rings 
[14], Pure haeckelite structures and their hybrids with graphene have a significantly 
lower thermal conductivity than pure graphene [15]. Only few theoretical studies have 
investigated electronic transport in selected graphitic carbon nanostructures with non- 
hexagonal rings including hybrid graphene-haeckelite structures [16, 17] and haeckelite 
nanotubes [18, 19, 20] 

Here we combine ab initio density functional theory (DFT) structural studies 
with DFT-based nonequilibrium Green function calculations to investigate how the 
presence of non-hexagonal rings affects electronic transport in graphitic structures. We 
find that infinite monolayers, finite-width nanoribbons and nanotubes formed of 5-8 
haeckelite with only 5- and 8-membered rings are generally more conductive than their 
graphene-based counterparts. Presence of haeckelite defect lines in the perfect graphitic 
structure, a model of grain boundaries in CVD-grown graphene, increases the electronic 
conductivity and renders it highly anisotropic. 

Haeckelites [8, 9, 10, 11, 12, 13] consist of periodic 2D arrangements of non- 
hexagonal rings of sp 2 bonded carbon atoms. Even though these structures have 
not been synthesized yet on a large scale, similar atomic arrangements have been 
observed (i) in 5-7 and 5-8 defect lines forming the in-plane interface between 
adjacent graphene flakes [14, 21, 22], (ii) in a vitreous atomic network formed during 
electron-beam irradiation of graphene [23], and (in) in 5-7 ring structures filling 
graphene nanoholes during the healing process of these defects [24], Most theoretical 
studies have focused on the equilibrium structure, stability and growth stability of 
haeckelites [8, 9, 10, 11, 12, 13, 17, 25, 26, 27, 28, 29] and found these systems to 
be either metallic, semi-metallic or semiconducting [12, 13, 27]. 

With continuing interest in CVD-grown graphene as an electronic material, 
increased attention must be given to carrier scattering at haeckelite-like grain boundaries 
connecting defect-free graphene regions. The most plausible model geometry to 
investigate charge transport in polycrystalline graphene is that of interconnected strips 
of haeckelite and graphene. A consistent picture should be obtained by comparing 
the effect of different structural arrangements and widths of haeckelite and graphene 
strips on the conductance and its anisotropy. As a counterpart to transport studies in 


Electronic Structure and Transport in Graphene/Haeckelite Hybrids: An Ab Initio Study 3 

graphene nanoribbons and nanotubes, we present corresponding results for haeckelite 
nanoribbons and nanotubes. 

2. Methods 

To gain insight into the equilibrium structure, stability and electronic properties of 
haeckelite structures, we performed DFT calculations as implemented in the SIESTA 
code [30]. Infinite 2D layers and ID ribbons and nanotubes were separated by 10 A 
thick vacuum regions in a 3D periodic arrangement. We used the Ceperley-Alder [31] 
exchange-correlation functional as parameterized by Perdew and Zunger [32], norm- 
conserving Troullier-Martins pseudopotentials [33], and a double-^ basis including 
polarization orbitals. The reciprocal space was sampled by a fine grid [34] of at least 
8x8x1 /.’-points in the Brillouin zone of the 2D primitive unit cell and its equivalent 
for ID structures or larger 2D supercells. We used a mesh cutoff energy of 180 Ry 
to determine the self-consistent charge density, which provided us with a precision in 
total energy of <2 meV/atom. All geometries have been optimized using the conjugate 
gradient method [35], until none of the residual Hellmann-Feynman forces exceeded 
10 -2 eV/A. 

Electronic transport properties were investigated using the nonequilibrium Green’s 
function (NEGF) approach as implemented in the TRAN-SIESTA code [36]. Ballistic 
transport calculations for optimized structures were performed using a single-^ basis 
with polarization orbitals, a 180 Ry mesh cutoff energy, and the same /.-point grid [34] 
as used for structure optimization. 

3. Results and discussion 

3.1. 5-8 haeckelite 

A perfect 2D monolayer of 5-8 haeckelite containing only 5- and 8-membered rings is 
shown in Fig. 1(a). The optimum rectangular unit cell is spanned by the Bravais lattice 
vectors ai=4.87 A in the ^-direction and a 2 =6.93 A in the ^-direction. The 5-8 haeckelite 
structure is about 0.36 eV/atom less stable than graphene, which is comparable to the 
stability of narrow carbon nanotubes. Due to this relatively high stability, we expect 
5-8 haeckelite structures to coexists with graphene at grain boundaries. 

The electronic band structure of 5-8 haeckelite is presented in Fig. 1(b). In contrast 
to semimetallic graphene, 5-8 haeckelite is metallic and has a finite electronic density of 
states at the Fermi level. Whereas the Fermi surface of graphene consists of 6 isolated 
/[■-points, that of 5-8 haeckelite is a line of finite length that intersects the T-Y high- 
symmetry line, as seen in Fig. 1(b). 

The results of our quantum transport calculation for this system are shown in 
Fig. 1(c). Besides the improved conductivity over graphene, suggested by the increased 
density of states at the Fermi level, we find the conductivity to be also anisotropic, as 
expected when considering the atomic arrangement in Fig. 1(a). We find the electrical 
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Figure 1 . (Color online) (a) Atomic structure and (b) electronic band structure of 5-8 
haeckelite with only 5- and 8-membered rings. The primitive unit cell is indicated by 
the shaded region, (c) Electronic conductance G jGo along the horizontal ^-direction, 
shown by the solid red line, and the vertical //-direction, shown by the dashed green 
line. G is normalized by the width of the unit cell normal to the transport direction. 
E = 0 corresponds to carrier injection at Ep- 


conductance along the y direction to be much higher than along the .'//-direction, and 
even observe a very narrow transport gap at Ep. These findings are consistent with 
a very anisotropic Fermi surface that crosses the T-Y, but not the F-X high-symmetry 
line in Fig. 1(b). 

3.2. 5-8 haeckelite nanoribbons and nanotubes 

Finite-width graphene nanoribbons and carbon nanotubes have received wide attention, 
since - unlike infinite graphene monolayers - some of these structures display sizeable 
band gaps. In analogy to these structures, we also studied quantum transport in ID 
haeckelite nanoribbons (h-NRs) and nanotubes(h-NTs). In Fig. 2, we present our results 
for h-NRs with different widths W that are passivated by hydrogen at the edge. The 
atomic structure of the three narrowest haeckelite nanoribbons are shown in Figs. 2(a-c). 
We find the optimum lattice constant T~4.9 A to be nearly independent of the width 
W. We also note that structures in Fig. 2(a) and 2(c) have mirror symmetry, whereas 
that in Fig. 2(b) does not. 

As seen in panels (d-i) of Fig. 2, we found the electronic structure and conductance 
of h-NRs to depend sensitively on the ribbon width. Structures with an odd number 
of 8-membered rings in the unit cell, depicted in Figs. 2(a) and 2(c), are narrow-gap 
semiconductors, with the fundamental band gap decreasing with increasing width from 
E g = 0.05 eV in Fig. 2(a,d) to E g = 0.02 eV in Fig. 2(c,f). Structures with an even 
number of 8-membered rings in the unit cell are all metallic. One example with two 
8-membered rings in the unit cell is shown in Fig. 2(b,e). Quantum conductance G of 
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Figure 2. (Color online) Atomic structure, electronic band structure and transport 
properties of haeckelite nanoribbons passivated by hydrogen at the edge. Atomic 
structure of haeckelite nanoribbons with various width are shown in (a), (b), (c). The 
corresponding electronic band structure and density of states is shown in (d), (e), (f), 
and ballistic transport conductance G/Go in (g), (h), (i). The density of state is given 
per atom. 


the three h-NRs in units of the conductance quantum Go is shown in Fig. 2(g-i). Of 
most interest is the presence or absence of a transport gap at E — 0, corresponding 
to carrier injection at E F . As expected, the semiconducting nanoribbons depicted in 
Figs. 2(a) and 2(c), have also a finite transport gap, seen in Figs. 2(g) and 2(i). The 
metallic nanoribbon in Fig. 2(b) does not have a transport gap at E = 0 according to 
Fig. 2(h). 

Similar to graphene-based carbon nanotubes, we can construct 5-8 haeckelite-based 
nanotubes, and characterize them by the chiral index ( n , m) in analogy to carbon 
nanotubes. We present the structure and electronic properties of two representative 
5-8 haeckelite nanotubes, h-NT(0,2) and h-NT(2,0), in Fig. 3. Both the side and the 
end-on view of these nanotubes in Fig. 3(a,b) indicates that their optimum cross-section 
is not as round and their surface not as smooth as that of their graphitic counterparts, 
owing to the presence of 5- and 8-membered rings. We found the ultra-narrow h- 
NT(0,2) and h-NT(2,0) nanotubes to be stable, but highly strained. The stability of 
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Figure 3. (Color online) Structural and electronic properties of haeckelite nanotubes. 
Atomic structure of (a) h-NT(0,2) and (b) h-NT(2,0) nanotubes. The length of the 
unit cells is denoted by L and the width by W. The corresponding electronic band 
structure and density of states is shown in (c), (d), and ballistic transport conductance 
G/Gq in (e), (f). The density of state is given per atom. 


the narrower h-NT(2,0) is lower by 0.81 eV/atom and that of the wider h-NT(0,2) 
is lower by 0.45 eV/atom with respect to the planar haeckelite structure depicted in 
Fig. 1(a). 

The electronic band structure and density of states of the h-NT(0,2) nanotube, 
shown in Fig. 3(c), indicates that h-NT(0,2) is metallic and has a non-zero density of 
states at the Fermi energy. As expected, also the calculated quantum conductance, 
shown in Fig. 3(e), indicates non-vanishing quantum conductance at E — 0. Quite 
different are the electronic structure and quantum conductance results for the h- 
NT(2,0) nanotube, shown in Fig. 3(d) and 3(f), which display a very small fundamental 
and transport gap of < 0.1 eV. Unlike in graphene-based carbon nanotubes, there 
is no general rule to predict whether a given h-NT nanotube should be metallic or 
semiconducting. 
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Figure 4. (Color online) Structure and electronic properties of hybrid 
haeckelite/graphene structures h-(n,m) consisting of graphene strips interconnected 
by 5-8 haeckelite strips, with n refereeing to the width of graphene and m to that of 
haeckelite. The upper panels depict the atomic structure of (a) h-(1,1), (b) h-(3,l) 
and (c) h-(5,l) hybrid structures with m = 1. The corresponding electronic densities 
of states are shown in the lower panels (d), (e) and (f). The length of the unit cells is 
denoted by L and the width by W. The density of state is given per atom. 


3.3. Hybrid haeckelite-graphene structures 

As a model of haeckelite-like grain boundaries connecting graphene grains in 
polycrystalline graphene samples, we construct hybrid haeckelite-graphene structures 
consisting of strips of 5-8 h-NRs of various width inter-connecting bare zigzag graphene 
nanoribbons with different widths. The hybrid systems, identified as h-(n, rn), are 
characterized by the number m of 8-membered rings per unit cell and the number n of 
hexagonal rings across the width of the unit cell. We explored two families of hybrid 
structures with m = 1 and m = 2, but different values of n. These studies let us explore 
how the density of 5-8 line defects should affect the electronic and transport properties 
of the hybrid structures. 

The optimum atomic arrangement and electronic structure of h-(n, 1) hybrids, with 
n — 1,3, 5, is shown in Fig. 4. We found all the structures to be metallic, as indicated by 
the peak in electronic densities of states at the Fermi energy in Figs. 4(d-f). The peak 
is associated with the 5-8 defects and decreases in strength with decreasing fraction of 
these defects in the structure, as the graphene strips become wider. 
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Figure 5. (Color online) Structure and electronic properties of hybrid 
haeckelite/graphene structures h-(n,m) consisting of graphene strips interconnected 
by 5-8 haeckelite strips, with n refereeing to the width of graphene and m to that of 
haeckelite. The upper panels depict the atomic structure of (a) h-(2,2), (b) h-(4,2) 
and (c) h-(6,2) hybrid structures with m = 2. The corresponding electronic densities 
of states are shown in the lower panels (d), (e) and (f). The length of the unit cells is 
denoted by L and the width by W. The density of state is given per atom. 


The corresponding results for h-(n, 2) hybrids, with wider haeckelite strips (m = 2) 
and also wider graphene strips characterized by n = 2,4, 6, are shown in Fig. 5. Similar 
to the m = 1 family presented in Fig. 4, also the h-(n, 2) structures are all metallic, 
but the peak in the density of states peak near Ep does not appear as sharp as in the 
h-(n, 1) structures. Similar to h-(n, 1) hybrids, the effect of the defect line diminished 
with increasing n. 

For the sake of illustration, we included the density of states of pristine graphene in 
the same energy range and observe that presence of 5-8 haeckelite line defects increases 
the density of states near the Fermi level, benefitting conductivity. 

We note that the narrow haeckelite strips within the h-(n, 1) hybrid structures 
correspond to the free-standing, but hydrogen-passivated nanoribbon presented in 
Fig. 2(a). Similarly, also the wider haeckelite strips within the h-(n, 2) hybrid structures 
have a counterpart in the free-standing, but hydrogen-passivated nanoribbon presented 
in Fig. 2(b). Comparison between the density of states of the free-standing haeckelite 
nanoribbons in Fig. 2(d,e) and the densities of states in Figs. 4(d-f) and Figs. 5(d-f) 
indicates that the electronic structure of free-standing and embedded nanoribbons is 
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Figure 6. (Color online) Electronic conductance G of hybrid haeckelite/graphene 
structures h-(n,m) consisting of graphene strips interconnected by 5-8 haeckelite strips. 
Results are presented for (a) h-(l,l), (b) h-(3,l), (c) h-(5,l), (d) h-(2,2), (e) h-(4,2), (f) 
h-(6,2) with different widths n of the graphene strips and m of the haeckelite strips. 
The structures and transport directions are defined in Figs. 4(a-c) and 5(a-c). G is 
normalized by the width of the unit cell normal to the transport direction. E = 0 
corresponds to carrier injection at Ep. 


very different. 

The increasing similarity of h-(n, m) structures with increasing width n of the 
graphene strips and defect-free graphene is also reflected in the quantum conductance 
of these systems, shown in Fig. 6. Unlike the density of states, quantum conductance 
is anisotropic. Especially for large values of n, transport along the defect lines will 
increasingly resemble that of graphene superposed with that of the conducting defect 
lines acting as quantum conductors or metal wires. We conclude that lines of non- 
hexagonal rings at grain boundaries in polycrystalline graphene should enhance the 
conductance of this system over pristine semimetallic graphene. 

4. Conclusions 

In conclusion, we have combined ab initio density functional theory (DFT) structural 
studies with DFT-based nonequilibrium Green function calculations to study how the 
presence of non-hexagonal rings affects electronic transport in graphitic structures found 
in polycrystalline graphene. We found that infinite monolayers, finite-width nanoribbons 
and nanotubes formed of 5-8 haeckelite with only 5- and 8-membered rings are generally 
more conductive than their graphene-based counterparts. Presence of haeckelite defect 
lines in the perfect graphitic structure, a model of grain boundaries in CVD-grown 
graphene, increases the electronic conductivity and renders it highly anisotropic. 
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